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Abstract. We provide details on a recent solution of the nucleon's covariant Faddeev equation in an explicit 
three-quark approach. The full Poincare-covariant structure of the three-quark amplitude is implemented through 
an orthogonal basis obtained from a partial-wave decomposition. We employ a rainbow-ladder gluon exchange 
kernel which allows for a comparison with meson Bethe-Salpeter and baryon quark-diquark studies. We describe 
the construction of the three-quark amplitude in full detail and compare it to a notation widespread in recent 
publications. Finally, we discuss first numerical results for the nucleon's amplitude. 



1 Introduction 

Abundant information on the structure of the nucleon is 
available from experiments, and here especially from elec- 
troweak probes at all energy scales. It is presently both an 
issue of paramount importance and a formidable theoret- 
ical challenge in contemporary hadron physics to under- 
stand the nucleon and its structure in terms of quarks and 
gluons which are the elementary degrees of freedom of 
quantum chromodynamics (QCD), i.e. the quantum field 
theory of the strong interaction. 

The formalism which has been developed to treat the 
three-body bound-state problem has a longstanding history 
which dates back to the original work by Faddeev Q. 
Non-relativistic Faddeev equations have found widespread 
application in the description of three-nucleon systems, 
see Ref. ||2) for an overview. The covariant generalization 
of the Faddeev equation to the three-body analogue 
of a two-body Bethe-Salpeter equation (BSE) Q was 
formulated in Refs. B4I5II : a comprehensive introduction 
can be found in [6 |. Within such a framework, the equation 
describes the baryon as a bound state of three spin- 72 
valence quarks where the interaction kernel comprises 
two- and three-quark contributions. 

The formalism of QCD's Dyson-Schwinger equations 
(for recent reviews, see e.g. [7 8]) provides a way to em- 
bed the covariant three-quark Faddeev equation in a con- 
sistent quantum-field theoretical setup. The dynamical in- 
gredients in the equation can then be treated in perfect cor- 
respondence with studies of quark and meson properties 
as well as related aspects of QCD. A solution of the equa- 
tion relies upon knowledge of the dressed quark propaga- 
tor and the three-quark kernel as well as a specification of 
the Poincare-covariant baryon amplitude. The relativistic 
spin structure of the latter has been explored in [9 lOj and 
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described in the light-front formalism in H 1 1 1 1 21 1311 41 . A 
complete classification according to the Lorentz group and 
the permutation group §3 was derived in |fl5l in terms of 
covariant three-spinors. 

The biggest obstacle on the way to a direct numer- 
ical solution of the three-body bound-state equation is 
its complexity. Upon implementing perturbative quark 
propagators it has been studied, for instance, in the 
works of Refs. 1161171 . in the context of a three-body 
spectator approximation [18|, or a Salpeter equation with 
instantaneous forces (5J. The corresponding equation 
of a scalar three-particle system with scalar two-body 
exchange based on the Wick-Cutkosky model [ 19 20] 
was recently investigated and compared to the light-front 
approach lt2T1 . 

A different kind of simplification can be achieved 
by considering diquark correlations in the nucleon, 
see e.g. Il22l for an overview. While maintaining full 
Poincare covariance, the quark-diquark model traces the 
nucleon's binding to the formation of colored scalar- and 
axialvector diquarks, thereby simplifying the Faddeev 
equation to a quark-diquark BSE. This strategy has 
been applied to investigate nucleon and A properties 
H23I24I25I26I27I28I29I30I3T1 . For more de tailed s tudies of 
diquark properties in this approach, see [32 33 1. Further 
support for the diquark concept has been provided by a 
study of diquark confinement in Coulomb-gauge QCD 
G4). 

In a recent study we have reported the first fully 
Poincare-covariant computation of the nucleon's mass and 
Faddeev amplitude beyond the quark-diquark approxima- 
tion [ 35 1 . The numerical solution of the Faddeev equation 
is performed after truncating the interaction kernel to a 
dressed gluon-ladder exchange between any two quarks, 
thereby enabling a direct comparison with corresponding 
meson studies as well as earlier investigations of baryons 
in the quark-diquark model. In the present work we pro- 
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vide details of that calculation with regard to the covariant 
structure of the nucleon amplitude and its decomposition 
into Dirac tensors, as well as details on the numerical 
solution of the Faddeev equation. 



2 Faddeev amplitude and equation 

In QCD baryons appear as poles in the three-quark scat- 
tering matrix. This allows one to derive a relativistic three- 
body bound-state equation: 



and use the value f = 1/3 in connection with the momen- 
tum p; this value maximizes the upper boundary for the 
nucleon mass with respect to restrictions arising from the 
quark propagator's singularity structure ll36l . The quark 
propagators S depend on the internal quark momenta kj = 
Pi - k and kj = pj + k. The internal relative momenta are 
given by: 
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where W is the bound-state amplitude defined on the 
baryon mass shell. The three-body kernel Kq) comprises 
a three-quark irreducible contribution and the sum of 
permuted two-quark kernels whose quark-antiquark 
analogues appear in a meson BSE. The subscript a denotes 
the respective accompanying spectator quark. 

The observation of a strong attraction in the S U(3)c 
antitriplet qq channel has been the guiding idea for the 
quark-diquark model, namely that quark-quark corre- 
lations provide important binding structure in baryons. 
This motivates the omission of the three-body irreducible 
contribution from the full three-quark kernel. The resulting 
covariant Faddeev equation includes a sum of permuted 
qq kernels (cf. Fig.Q~|l: 

?W(P, q,P) = J] f K ( :i myY ?W*(p (a) , q (a \ P) , 

(2) 

where denotes the renormalization-group invariant 
products of a qq kernel and two dressed quark propagators: 
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Here, {a, b, c] is an even permutation of {1,2, 3} and linked 
to the respective Dirac index pairs. 

The spin-momentum part of the full Poincare-covariant 
nucleon amplitude ¥ a p y s(p,q,P) is a spin-72 four-point 
function with positive parity and positive energy: it car- 
ries three spinor indices {a,/3, y) for the involved valence 
quarks and one index 6 for the spin- 1/2 nucleon. The am- 
plitude depends on three quark momenta p\, p2, pi which 
can be reexpressed in terms of the total momentum P and 
two relative Jacobi momenta p and q, where P 2 = -M 2 is 
fixed. They are related via (cf. Fig.Q~|i: 



p = (1 - 0P3-(Pd, 
P2~P\ 

P = Pi + P2 + P3 , 



P\ = —Q P , 
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where we abbreviated pd :- p\+ p2- We have chosen equal 
momentum partitioning 1 /2 for the relative momentum q 



The nucleon amplitude can be decomposed into 64 
Dirac structures: 



64 



VapysiP, q> P ) = 2_j fk(p 2 , q 2 , fel) T k(P, <7, P)afly6 , (6) 



k=\ 



where the amplitude dressing functions depend on the 
five Lorentz-invariant combinations 

p 2 , q 2 , Zo=pf-q~T, Zi=p'P, Z2 = q ■ P ■ (7) 

Here, a hat denotes a normalized 4-vector and pj = Tp V p v 
a transverse projection with respect to any four momentum 
P: Tp" = 5^ - Pfp*. We abbreviated the angular variables 
by the shorthand notation {z} = {zo, Z\ , Zi}- The Dirac struc- 
tures Tk(p, q, P) will be explained in SectionH 



3 Rainbow-ladder truncation 

To proceed with the numerical solution of the Faddeev 
equation, we need to specify the quark-quark kernel and 
the dressed quark propagator S (p) which appear in Eq. (0. 
They are related via the axial-vector Ward-Takahashi iden- 
tity (AXWTI) which encodes the properties of chiral sym- 
metry in connection with QCD. Its satisfaction by the in- 
teraction kernels in related equations guarantees the cor- 
rect implementation of chiral symmetry and its dynamical 
breaking, leading e.g. to a generalized Gell-Mann-Oakes- 
Renner relation valid for all pseudoscalar mesons and all 
current-quark masses 1 37 38 1. In particular the pion, being 
the Goldstone boson related to dynamical chiral symmetry 
breaking, becomes massless in the chiral limit, indepen- 
dent of the details of the interaction. Specifically, we de- 
scribe the qq kernel by a dressed gluon-ladder exchange: 



47Ta(k 2 ) 



(8) 



where k is the gluon momentum and Z2 the quark renor- 
malization constant. By virtue of the AXWTI, the kernel 
must also appear in the quark Dyson-Schwinger equation 
whose solution defines the renormalized dressed quark 
propagator: 

S~l{p) = Z 2 {ip + m) af) + f % aaW {k)S a ,p,{q) , (9) 

Jo 
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Fig. 1. Faddeev equation ifJJ in rainbow-ladder truncation. 



where the bare quark mass m enters as an input and k = 
p - q. The inherent color structure of the kernel leads to 
prefactors 2/3 and 4/3 for the integrals in Eqs. (O and (O, 
respectively. 

Eqs. (J8H9]) define the rainbow-ladder (RL) truncation 
which has been extensively used in DSE studies of mesons 
and baryons in the quark-diquark model, e. g. [39 27j and 
references therein. The non-perturbative dressing of the 
gluon propagator and the quark-gluon vertex are absorbed 
into an effective coupling a(k 2 ) for which we adopt the 
ansatz 1401411 

a(k 2 ) = 7T7] 7 J^J e~^) + auv(k 2 ) . (10) 



The second term 

auv(k 2 ) = 



In ^e 2 -\+{\+k 2 IA 2 QCD f 



(11) 



where y m = 12/(1 INc - 2JV/) is the anomalous dimension 
of the quark propagator, reproduces the logarithmic 
decrease of QCD's one-loop perturbative running cou- 
pling and vanishes at k 2 = 0. In our calculation we use 
Aqcd = 0.234 GeV and y m = 12/25 which corresponds to 
Nf = 4. The infrared contribution (the first term in ( TTOl i) is 
parametrized by an infrared scale A and a dimensionless 
parameter 77 and yields the non-perturbative enhancement 
at small and intermediate gluon momenta necessary to 
generate dynamical chiral symmetry breaking and hence 
a constituent-quark mass scale. ({A, rf\ and the infrared 
parameters used in [41] are related by C — (A/A,) 3 and 
cj = rf l AIA,, with A, = 1 GeV.) 

Beyond the present truncation, corrections arise from 
pseudoscalar meson-cloud contributions which provide a 
substantial attractive contribution to the 'quark core' of 
dynamically generated hadron observables in the chiral 
regime and vanish with increasing current-quark mass, but 
also from non-resonant contributions due to the infrared 
structure of the quark-gluon vertex. To anticipate such 
corrections we exploit the freedom in adjusting the input 
scale A. We adopt two different choices established in the 
literature in the context of n and p properties [41|: 

Setup A is determined by a fixed scale A = 0.72 GeV, 
chosen in [40 1 to reproduce the experimental pion decay 
constant and the phenomenological quark condensate. Cor- 
responding results are therefore aimed in principle at a 
comparison to experimental data for meson and baryon 
properties (see H39I29H and references therein). Setup B de- 
fines a current-mass dependent scale which is deliberately 



inflated close to the chiral limit, where A ~ 1 GeV [41 1. 
It is meant to describe a hadronic quark core which must 
subsequently be dressed by pion-cloud effects and other 
corrections. As a result, ji, p, N and A observables are con- 
sistently overestimated, but (with the exception of the A- 
baryon) compatible with quark-core estimates from quark 
models and chiral perturbation theory (for a detailed dis- 
cussion, see [41 28 29|). Irrespective of the choice of A, 
hadronic ground-state properties have turned out to be in- 
sensitive to the value of r] in a certain range [40 39). Con- 
sequently, with Eqs. (TSl-FTOb and A, the input of the Fad- 
deev equation is completely specified with all parameters 
already fixed to meson properties. 



4 Structure of the nucleon amplitude 

With the kernel K of the Faddeev equation (f2]) determined, 
we still need to find expressions for the basis elements 
which constitute the nucleon amplitude according to 
Eq. ©. A general spinor four-point function which 
depends on 3 independent momenta involves 128 com- 
ponents of positive parity. These can be conveniently 
expressed through tensor products of two Dirac matrices, 
for which we adopt the following notation: 

(A ® B) a fr yS = A a pB y s , 
(Ai ® Bi)(A 2 ® B 2 ) = (AiA 2 ) ® (BiB 2 ) , 
(A ® B) T =A T ®B T , 
Tr {A ® B] = Tr{A)Tr{S}. 

We define a charge-conjugated four-point function via 



{A®B){p,q,P) := 

= {C® C)(A ® B) T (-p, -q, -P)(C ® Cf = (12) 

= A(p, q, P) ® B(p, q, P) , 

where C = y 4 y 2 is the charge-conjugation matrix. A 
further classification into the 64-dimensional subspaces 
of positive and negative energy, describing the nucleon's 
(l/2) + and (l/2) + states, is obtained by attaching the 
respective projectors A ± (P) = (1 + f)/2 to the nucleon leg 
in the tensor product: A ® BA ± . The negative-parity states 
(l/2)~ and (l/2)~ constitute another 128-dimensional 
vector space generated from basis elements of the type 
Ay 5 ®BA ± . 

To construct a suitable basis for the positive-parity and 
positive-energy nucleon, we start from the expressions 



if(P) = A r (P)y 5 C®A + (P), 



(13) 
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with r - + according to the projectors A r (P). In the quark- 
diquark model, Q + + Q~ — 75 C ® A + is related to the dom- 
inant scalar-diquark part in the nucleon amplitude. To in- 
clude the momentum dependence of the amplitude, it is 
convenient to choose a set of momenta \p~r, q~,, P] which 
are orthonormal with respect to the Euclidean metric, i.e. 



p T =q t = P = 1, 
Pt -q~t =P~t ■ P = <i~f P 
This can be realized via 



0. 
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Consider now the four quantities 

r i ( P ,q,p)e{t, \ (hJ t ], h, 4,} 



(14) 
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(16) 



where, because of Eq. (fl4l . one has 4 [pr,4t] = Pt it ■ The 
linear span of the 8 matrices ri{p,q,P) A'\P) equals the 
linear span of the basis elements 

whereas /",(/?, q, P) A + (P) reduces this set to 

{t,p,4,P4)A + (18) 

since P A + = A + . Hence the quantities (r, ® rj) Q' define 
32 linearly independent basis tensors. 

A complete orthogonal basis for the nucleon amplitude 
is given by the 64-dimensional set {Sr., Pj\} defined by 



'%<J>,q,P)). 
P r {J (p,q,P) 



(19) 



with i, j = 1 ... 4. The symbols S and P were chosen to 
reflect a combination of two scalar or pseudoscalar covari- 
ants whose product again exhibits positive parity. All pos- 
sible further basis elements, e.g. 



'^j(p,q,P)) 
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(r i ®r j )Q r , (20) 


UljiP^P), 






$0^' , 





where jj = Tp V y v and = -i[y l T ,y T ]/2, linearly de- 
pend on those in Eq. ( fT9] i; respective relations are given 
in Table Q] For instance, the dominant contributions to the 
Faddeev amplitude are the elements 

r=± 

f T C ® f T y 5 A + = £ = J] (r S^ 2 + P 33 + PJ4) . 

r=± r=± 

In the quark-diquark model, these correspond to scalar- 
scalar and axialvector-axialvector combinations of diquark 
and quark-diquark amplitudes for either of the three dia- 
grams appearing in the Faddeev equation in Fig.Q] 



A partial-wave decomposition (see Tableland App.lBl) 
leads to linear combinations of the {S^, P[.} as eigenstates 

of quark-spin and orbital angular momentum operators S 2 
and L 2 in the nucleon rest frame. The 64 basis covariants 
(32 each for total quark spin 5=1/2 and s = 3/2, respec- 
tively) can be arranged into sets of 8 s-waves (I = 0), 36 
p-waves {I — 1), and 20 d- waves (/ = 2) which we denote 
collectively by 



1 ® 1 
y 5 ® y 5 



T u Q r , 



(21) 



with i, j = 1 ... 4 and r — +. The 7y carry the relative- 
momentum dependence and are defined in Table [2] For in- 
stance, the eigenstates of s = l /2 correspond to i = 1,2 and 
read 



1 



1 Pi 



(22) 



whereas the s = 3 /2 eigenstates (z = 3,4) are more com- 
plicated. The relations between the X[-j and the basis ele- 
ments {Sr P.j , A.j, \/.j} are stated in fable [2| as well. 

As illustrated in App.lCl the basis elements X 1 .. can be 
expressed in terms of quark three-spinors frequently used 
in the literature, e.g. Ref. [15|. Moreover, they satisfy the 
following orthogonality relation: 



; Tr {X; 7 X^,} = \ 0<$j)fiafy 0^f)c^,yS = S » 6 if 



(23) 



As an example, consider the scalar-scalar combinations: 
%(p, q,P)=C [Tii-p, -q, -P) A r (-P) 75C} T C T 

®c{rj(-p,-q,-P)A + (-p)} T c T = 

= (C T y 5 )(A r f i ) ®(A + fj), 

where we used the relations A r (P) = C A r (-P) T C T = 
A r (P) and Ti(p, q,P) = C r t (-p, -q, -P) T C T . As a result 
one obtains 



iTr{S^S^,} =±Tr{r i r i ,A r A r }Tr{r j r f A + } = 
=Su> S jf 5 rr > , 

since A r ' A r = 5 rf A r and Tr [T { r r A r ) = 2 8 a , . 

The Pauli principle requires the Faddeev amplitude 
to be antisymmetric under exchange of any two quarks. 
The Faddeev kernel is invariant under the permutation 
group §3. The eigenstates of the Faddeev kernel can hence 
be arranged into irreducible §3 multiplets 

V * 7 *'(l%)' (24) 

of which the first two (totally symmetric or antisymmetric) 
solutions are unphysical while the mixed-symmetry dou- 
blet constitutes the Dirac part of the nucleon amplitude. 
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Table 1. Relations between the basis elements Ay of Eq. J20I ) and (Sy, Py). The corresponding relations for Vy are obtained by in- 
terchanging Sy «-> Py, Similar dependencies hold for the Ty, e.g.: T u = -2A|,. The superscripts r = ± are not displayed for better 
readability. 
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Table 2. Orthonormal Dirac basis XJ jt of Eq. d21b constructed from a partial-wave decomposition. The first two columns denote the 
eigenvalues of total quark spin s and intrinsic orbital angular momentum / in the nucleon rest frame. The third and fourth columns 
define the relation between the X'. and the basis elements from Eq. l l 19120b . Each row involves 4 covariants; the superscripts r = ± are 
not displayed for better readability. The fifth column shows the momentum-dependent covariants Fy which appear in Eq. ( 121b : we have 
abbreviated pf — » p and q~, —* q for clarity. 
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Fig. 2. Zeroth Chebyshev moments of the dressing functions cor- 
responding to the dominant covariants in the Faddeev amplitude 
^V?,, plotted as a function of yp^ and yq^. The various curves 
represent the falloff in the variable -Jq^. 
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Fig. 3. First four Chebyshev moments in the variable z\ of the 
dressing functions associated with the amplitudes S,,, evaluated 
at q- = 0. 



Taking into account the flavor and color structure, the full 
Dirac-flavor-color amplitude reads 



¥(p,q,P) = [W M:n J M:n + WmJms}^ 



(25) 



where babc is the antisymmetric color-singlet wave func- 
tion and 

T M* = -^5 io~2 ® 1 . = --^= (rio-2 <8> cr (26) 

denote the isospin-1/2 flavor tensors involving the Pauli 
matrices <x,. A projection onto the proton or neutron flavor 
states involves a contraction of the rearmost flavor index 
with either of the two isospin basis states (1,0) or (0, 1). 

A flavor-dependent kernel in the Faddeev equation 
will mix fjviy, an d ^M s whose dominant contributions 
are given by and A| p respectively. However, since 
the rainbow-ladder kernel presently employed is flavor- 
independent and we consider only equal quark masses, 
the equations for the Dirac amplitudes W\\^ and ¥m s m 
Eq. d25l l decouple because of the orthogonality of the two 
flavor tensors 1 m a an d ~^M S '- 



V = K ) V 



~ ^(3) ¥m m , 
'PjHs - ^(3) ^ Ms- 



ill) 



Hence one obtains two degenerate solutions of the Faddeev 
equation, where by virtue of the iterative solution method 
the symmetry of the start function determines the symme- 
try of the resulting amplitude. 

Similarly to the analogous case of a diquark ampli- 
tude, the quark exchange symmetry does not reduce the 
number of Dirac covariants since the dressing functions 
fkip 2 , q 2 , {z}) transform under the permutation group 
as well. For example, if we restrict ourselves to the 
eight momentum-independent covariants of Table [2] (with 
s - t/2, / = 0) and apply the §3 (anti-)symmetrizers, we can 
arrange these 8 basis elements into three mixed-symmetry 
doublets and one symmetric and one antisymmetric singlet 
which are shown in Table [3] If the coefficients fa were 
totally symmetric, e.g., by being constant or by depending 
only on certain symmetric combinations of p 2 , q 2 and {z} 
as derived in Ref. 03], the nucleon amplitude would be a 
linear combination of the six My\ and Ms basis elements 
in Tabled 



k 

AU 



AU 



k=] 



* 

Ai s ■ 



(28) 



k=\ 



Since the coefficients /* can appear in all symmetry repre- 
sentations, the inclusion of the remaining Dirac covariants 
ifffi and if/g is however necessary, and the same reasoning 
holds if all 64 basis elements are implemented. 



5 Results 

The explicit numerical implementation of the Faddeev 
equation is described in App. [A] The massive computa- 
tional demand in solving the equation primarily comes 



19 International IUPAP Conference on Few-Body Problems in Physics 



from the five Lorentz-invariant momentum combinations 
of Eq. (0 upon which the amplitudes depend. In analogy 
to the separability assumption of the nucleon amplitude in 
the quark-diquark model we omit the dependence on the 
angular variable Zo =pr ■ q~r but solve for all 64 dressing 
functions fk(p 2 , q 2 , 0,zu Zz)- 

The resulting nucleon masses at the physical pion mass 
in both setups A and B are shown in Table [4] As a conse- 
quence of Eq. $T7\ , the two states and Ym s emerge as 
independent solutions of the Faddeev equation. Both sep- 
arate equations produce approximately the same nucleon 
mass, where the deviation of ~ 2% is presumably a trun- 
cation artifact associated with the omission of the angle zo- 
For either solution typically only a small number of co- 
variants are relevant which are predominantly s-wave with 
a small p-wave admixture. The corresponding amplitudes 
for the mixed-antisymmetric solution are shown in Fig. [2] 
Comparing the relative strengths of the amplitudes allows 
to identify the dominant contributions: 

^Mn : ^{Sn, Vj , , Si 3 , V| 3 , X{J 3)1 }, 

A , , (29) 

r 

Fig. 3 displays the angular dependence in the variable 
Z\ through the first few Chebyshev moments of the ampli- 
tudes S^j which contribute to ^Mm- The angular depen- 
dence in the variable zi is small compared to z\. This is 
analogous to the quark-diquark model, where the depen- 
dence on the angle between the relative and total momen- 
tum of the two quarks in a diquark amplitude is weak. 

The evolution of and the p-meson mass from the 
BSE vs. m 2 is plotted in Fig. [4] and compared to lattice 
results. The findings for are qualitatively similar to 
those for m p : setup A, where the coupling strength is ad- 
justed to the experimental value of f n , agrees with the lat- 
tice data. This behavior can be understood in light of a re- 
cent study of corrections beyond RL truncation which sug- 
gests a near cancellation in the p-meson of pionic effects 
and non-resonant corrections from the quark-gluon vertex 
[42 1 . Setup B provides a description of a quark core which 
overestimates the experimental values while it approaches 
the lattice results at larger quark masses. 

A comparison to the consistently obtained quark- 
diquark model result exhibits a discrepancy of only ~ 5%. 
This surprising and reassuring result indicates that a 
description of the nucleon as a superposition of scalar 
and axial-vector diquark correlations that interact with the 
remaining quark provides a close approximation to the 
consistent three-quark nucleon amplitude. 



6 Conclusions and outlook 

We have provided details on a fully Poincare-covariant 
three-quark solution of the nucleon's Faddeev equation. 
The nucleon amplitude which is generated by a gluon 
ladder-exchange is predominantly described by s- and 
p-wave Dirac structures, and the flavor independence 



Table 4. (adapted from Ref. [35]) Nucleon masses obtained from 
the Faddeev equation in setups A and B and compared to the 
quark-diquark result. The rj dependence is indicated for setup B 
in parentheses. 





Q-DQ|29| 


Faddeev 


Faddeev (Ms) 


Setup A 


0.94 


0.99 


0.97 


Setup B 


1.26(2) 


1.33(2) 


1.31(2) 




0.0 0.1 0.2 0.3 0.4 0.5 
m n 2 [GeV 2 ] 

Fig. 4. (adapted from Ref. 1 35 1) Evolution with m 2 n of m p and M N 
compared to lattice data ( 143 1441 for M N and 145 1461 for m p ). The 
quark-diquark model result for M N 1 28 29 1 is plotted for com- 
parison. Dashed and dashed-dotted lines correspond to setup A; 
the solid line for m p and the bands for M N (mixed-antisymmetric 
solution) are the results of setup B, where the variation with rj is 
explicitly taken into account. Dots denote the experimental val- 
ues. 



of the kernel leads to a mass degeneracy. The resulting 
nucleon mass is close to the quark-diquark model result 
which stresses the reliability of previous quark-diquark 
studies. 

Due to the considerable computational efforts involved, 
more results and an in-depth investigation with regard to 
the complete set of invariant variables will be presented 
in the future. Further extensions of the present work will 
include an analogous investigation of the zl-baryon, more 
sophisticated interaction kernels, e.g. in view of pionic cor- 
rections, and ultimately a comprehensive study of baryon 
resonances. 
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A Numerical implementation 

Similar to the analogous case of a two-body Bethe-Salpeter 
equation, the Faddeev equation (f2]i can be viewed as an 
eigenvalue problem for the kernel K^y. 

K (3) (P 2 ) Wi = Ai(P 2 ) Wi , (30) 

where P is the total momentum of the three-quark bound 
state and enters the equation as an external parame- 
ter. Upon projection onto given quantum numbers, the 
eigenvalues of K@y constitute the trajectories Aj(P 2 ). An 
intersection Ai(P 2 ) — 1 at a certain value P 2 = -M 2 
reproduces Eq. (f2]l and therefore corresponds to a potential 
physical state with mass M,. The largest eigenvalue Aq 
represents the ground state of the quantum numbers under 
consideration and the remaining ones Ai>\ its excitations; 
the associated eigenvectors ¥-, are the bound-state am- 
plitudes. (Note that in this context one has to keep in 
mind the possibility of anomalous states in the excitation 
spectra of bound-state equation solutions flTl .) To obtain 
the ground-state solution, Eq. ( f30b is solved via iteration 
within a 'guess range' P 2 e {-M 2 ^ m , -M^J, where M max 
is determined from the singularity structure of the quark 
propagator (see e.g. 136I48H ). Upon convergence of the 
eigenvalue Aq(P 2 ) the procedure is repeated for different 
P 2 until /l()(P 2 = -M 2 ) = 1, thereby defining the nucleon 
mass M. 

From a numerical point of view, it is advantageous to 
split the Faddeev equation for K@) = KS S into an equation 
for a 'wave function' <t> — SS ¥ which can be evaluated 
outside the loop integral, and a subsequent integration ¥ - 
K<t> that is carried out by calling the wave function <t> with 
the loop momenta as its arguments: 

& ia) (p,q,P) = S( Pb )S( Pc ) ¥{p,q,P), 

1 3 r 

¥(p, q, P) = g I K^\k) <P«(pM q (a \ P) . 

The index a - 1 ... 3 denotes the three permutations of 
the Faddeev kernel, and {a, b, c) is an even permutation of 
{1,2, 3}. For a complete Dirac basis the wave function 
can be projected onto the same basis elements as the am- 
plitude, cf. Eq. (|6]l, and we denote its dressing functions by 
ff' '. Exploiting the orthogonality relations ( 1231 yields cou- 
pled equations for the amplitude and wave function dress- 
ing functions: 




jt a \s) = Y J G<$(s)fj(s), 

i 



where the kernel K:^ and quark propagator matrix Gy are 

the matrix elements of K <a) and SS upon projection onto 
'outer' basis elements from the left and 'inner' basis ele- 
ments on the right. We abbreviated the five outer momen- 
tum variables by s = {p 2 ,q 2 ,Zo,Zi,Zi} and the respective 
inner variables by s (a) : they depend on s and the four com- 
binations t :- [k 2 , k ■ P, k ■ p, k ■ q}. 

Eq. (l3~TT l involves an iterated (multidimensional) 
matrix-vector multiplication, where in a straightforward 
implementation K and G would be computed in advance 
whereas the s (fl) dependence of the fj must be interpo- 
lated in each iteration step. The biggest obstacle in such an 
approach is the memory requirement of the kernel which, 
for reasonable accuracy, is of the order of Petabytes and 
hence far beyond the capacities of today's computing 
resources. In this respect it is advantageous to split the 
kernel K into the following three contributions: 

k^, = [Ki mr ys%\s, »] u , 02) 

where K; stems from the rainbow-ladder kernel (O 
modulo its Dirac structure and only depends on the 
gluon momentum k; K2 contains all Dirac traces but is 
momentum-independent; and K3 carries the remaining 
loop-momentum dependence inherent in the of Eq. ( [ToT l 
which enter the basis elements. Ki and K2 are independent 
of the baryon momentum P 2 and can be calculated and 
stored outside of the iterated matrix-vector multiplication. 
K3 is computed in each iteration step anew but includes 
only simple products of the loop momenta which, in the 
absence of any Dirac structure, are evaluated compara- 
tively quickly. A construction analogous to Eq. ( 1321 can 
be applied to the propagator matrix G'"' as well; however, 
here the situation is much less severe since no integration 
and hence no dependence on t is involved. 

Moreover, using the set of orthogonal momenta defined 
in Eq. dl~5b turns out to be extremely efficient: for instance, 
choosing the momentum alignment 

p= Jp^io, 0, Jl^z], zi), 

1 '_ (33) 

q = a/^{ ' V^ Zo ' Z2 } 

yields in the baryon's rest frame: 

?= {0,0, 0,1}, pi = {0, 0, 1, 0}, q, = {0, 1, 0, 0). 

As a consequence, all outer basis elements effectively do 
not depend on any momentum variable at all, and the mo- 
mentum dependence is solely carried by the inner basis el- 
ements which contribute to K^ fl) (.s, f). 

A combination of these strategies greatly reduces the 
memory demand to < 1 GB. The impact on the run time 
due to the on-the-fly computation of K3 is still slightly out- 
weighed by the time consumed to interpolate the dressing 
functions /* fl) inside the integral. To address this issue, we 
drop the dependence on the angular variable zo — pr ■ qj 
which, from analogy of the quark-diquark model, is ex- 
pected to be small. In addition we perform an expansion 
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into Chebyshev polynomials of the first kind in the re- 
maining angles z " , Zo that appear in the wave function 
coefficients inside the integral. The resulting run times are 
accessible by a parallel cluster. 

The different solutions I'm * an d ^Ms °f tne Faddeev 
equation are obtained by choosing suitable start functions 
for the amplitudes fi(s). To arrive at !Pai„, we start from the 
structure S|j and set all other amplitudes to zero; for Ym s 
we use the initial amplitude A|j . The full structure of the 
nucleon amplitude in either case is subsequently generated 
by iterating Eq. (13 It . 



and o^j = -i [jj, jj]/2. From these relations one can ver- 
ify that the XJ\ k are eigenfunctions of W 2 with j = 1/2 as 

wellas eigenfunctions of S 2 — > s(s+l)andL 2 — > /(7+1) 
which, in the nucleon's rest frame (where L 2 = L 2 and 
S 2 — S 2 ), assume the interpretation of total quark spin and 
intrinsic quark orbital angular momentum. Such a partial- 
wave decomposition yields 32 basis elements for each total 
quark spin 5=1/2 and s = 3/2, respectively, and a classi- 
fication into sets of 8 s waves (I = 0), 36 p waves (I — 1), 
and 20 d waves (I = 2). 



B Angular momentum decomposition 

The Dirac basis elements Xr.^ in Table|2]can be classified 
with respect to their quark-spin and orbital angular mo- 
mentum content in the nucleon's rest frame. Only the total 
angular momentum j - 1/2 of the nucleon is Poincare- 
invariant while the interpretation in terms of total quark 
spin and orbital angular momentum will differ in every 
frame. The spin is described by the Pauli-Lubanski opera- 
tor: 



2 



(34) 



where we chose the total momentum P to be normalized. 
J 1 " and W are the generators of the Poincare algebra satis- 
fying the usual commutation relations. The eigenvalues of 
the square of the Pauli-Lubanski operator, 



(35) 



define the spin j of the particle. For a system of three par- 
ticles with total momentum P and relative momenta p and 
q, the total angular momentum operator consists of the to- 
tal quark spin 5 and the relative orbital angular momentum 
L = L(p) + L(„). Subsuming them into Lorentz-covariant 
operators 

S" = \^ m Pp v {a° p ® 1 ® 1 + perm.) , 
L[ p) = Le» m Pp v (p a SP p - pPd*) 1 ® 1 ® 1, (36) 
L[ q) = L e » mfi P v (q a dP - <fd%) 1 ® 1 ® 1 

with W = + L'L + yields the identities 

5 2 = |l®l®l + i(cr^ v ®c7^ v ®l + perm.) , 
S -L {p) = { Y d v p (a^ ® 1 ® 1 + perm.) , 

L 2 p) = 2p -dp + yy^-fdl, 

L (pY L (q) = W ~ (P - q)(d P -d q ) , 
where we abbreviated 



p» = T p °'p v , cf p = T p !V d/dp Y , 



(37) 



C Multispinor notation 

The basis elements X r ^ k can be expressed in terms of the 
Dirac spinors U a (P), W(P) = y 5 t/ <r (P) which satisfy the 
free Dirac equation for a spin- 72 particle, i.e. which are 
eigenspinors of the positive and negative energy projectors 
A±. Normalized to U P (P) t/°"(P) = S pcr , they are given by 



U°~(P) = 



s + M 



2M 



w 

\ e+M 



(1,0) 
(0,1) 



(38) 



with e = Vf 2 + M 2 . A 64-dimensional basis for the nu- 
cleon wave function was constructed in Ref. [ 15 1 from the 
quark three-spinors UUU, VVU, VUV, UVV and their 
parity -reversed counterparts VVV, UUV, UVU, and VUU 
by equipping them with the 8 possible spin-up/down ar- 
rangements f |T, TTi, an d so on - Identities such as 

A + = t/ T t? T + U l U l , A + y 5 C = U 1 ^ - U^U 1 (39) 

allow to derive the following relations for the momentum- 
independent s-wave basis states of Table [2] 



-P- 



(U^U 1 - U l U^) t/ T . 
(yTyl _ ylyT) j/T ; 

t/T = (yTf/l _ yif/T)yT j 

(J/Tyl _ j/lyT) yT 



and 



A+ t/ T = (UW L + U [ W) t/ T - 2 UWW L , 

-A n t/ T = (v T v x + v x v T ) t/ T - 2 yiyW 1 , 

V+ t/T = (yTf/J. + ylj/T) yT _ 2 yTf/Tyl t 
-V- t/ T = (U^V 1 + U l Vl) V T - 2 C/TyTyl _ 

Due to the Poincare-covariant construction of Eqs. (119H21I I. 
the remaining 56 covariants depend on the relative mo- 
menta. For instance, the Dirac structure S n satisfies the 
relation 

;S B f/ T = (yTyi-yJ-yt) (([^ + i[p ^ h ) yi + [^] 3 yT) . 

(40) 

Using the special momentum alignment d33l leads to 



on 



(yTyl _ ylyT) yT 



(41) 
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and similar relations for the remaining basis elements in 
Table[2] From a three-spinor perspective, Eq. (HTh is at the 
same time the 'parity-flipped' counterpart of St v 

- (75 ® 75) (y T 5 ® t/ T ) = ( V T V x - V [ V r ) V T , (42) 

which would appear in the amplitude of a negative-parity 
nucleon (l/2~) and that of a positive-parity antinucleon 
(l/2 + ). Nevertheless this seemingly odd-parity structure 
still contributes to the (l/2 + ) state through Eq. d40l ), where 
the parity flip induced by the spinor replacement U — > V 
is saturated by an odd parity carried by the relative mo- 
mentum pr. As a consequence, indeed all 64 three-spinor 
combinations contribute to the nucleon's amplitude. 
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Table 3. Irreducible multiplets of the permutation group § 3 , constructed from the 8 covariants {S j l , P[ 1 , A\ l , VJ j }. 

^ = sr, ^ = 2 r Pn + s„ ^ = Zr(vj 1 -pj 1 )+2sr 1 ^ = E r (v; 1 + p; 1 )-2s ri 

^m s = A u ^M s =^ry r u -A- U ^, = 2 r r(V; i+ 3P; i ) + 2A u * s = £ r r(-\T u + 3PJ,) - 2 A~ 



